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I. INTRODUCTION 

Over 80 years ago, many physical systems were found 
to be governed by Bose-Einstein (BE) or Fermi-Dirac 
(FD) statistics [U El EJ SJ |5], as distinct from classi- 
cal degenerate Maxwell-Boltzmann (MB) statistics de- 
veloped in the 19th century [Bj [7J. The combinatorial 
weights of these statistics - the number of configurations 
(microstates) of each identifiable realization (complex- 
ion, macrostate) of the system - are widely given as [e.g. 
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where Stotai is the total dimensional (thermodynamic) 
entropy and k is the Boltzmann constant. Maximizing 
the entropy defined by |4]) ("MaxEnt"), subject to the 
constraints on the system, is therefore equivalent to iden- 
tifying the system with its most probable realization, a 
method of inductive reasoning which can be termed the 
maximum probability ("MaxProb") principle [51 [71 1131 1141 
IHailElinilTHllH]. Normally, @ is subject to the asymp- 
totic limit N —t oo (simplistically using the Stirling ap- 
proximation, lnm! sa mlnm — m, or more rigorously 
with Sanov's EO] theorem), to give the entropy function 
for each statistic [TlElElSH51EllZllHllSlIiniiniII21- 
Recently, Q has been applied directly to give the ex- 
act or non- asymptotic entropy function for each statistic 
[TBI [T71 [T51 [T9~] . enabling the application of statistical 
mechanics to systems containing small numbers of enti- 
ties. The non-asymptotic BE and FD entropy functions 
have peculiar information-theoretic properties, with im- 
portant implications for quantum mechanics |161 117j . 

The above weights were generalized by Brillouin [9, 
also in HUES to give: 



where i denotes each distinguishable level (outcome or 
state) of the system, from a total of s levels; m is the 
number of entities in each level i; N — Xa=i n i i s the 
total number of entities; and gi is the number of distin- 
guishable sublevels within each level i (the degeneracy or 
multiplicity) . In each case, the dimensionless entropy per 
entity H is obtained from the Boltzmann principle [3[7]: 
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where P denotes a correctly normalized probability dis- 
tribution and G = Yli=i 9* i s the total degeneracy. These 
are given in many references [e.g. 153] as the probabilities 

s^Nl/UUi^n-y 1 and (j^)" 1 respectively, i.e. 
without considering degenerate sublevels within each 
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level. In contrast, the following generalized distributions 
were given by Kapur |24j : 
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where (ft is the source probability 1 of an entity in level i. 
In the limits <7i<Cl and gi^>Ui, these imply: 
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(see later discussion). The Brillouin and Kapur sets of 
distributions are incompatible, and their range of validity 
and interrelationships have not been explored. 

Of great interest in recent years is the possible occur- 
rence of "anyon" particles, intermediate between bosons 
and fermions. Several intermediate statistics have been 
proposed. A long-standing philosophical approach is to 
consider the allocation of indistinguishable particles to 
distinguishable degenerate boxes, with a maximum of d 
particles per box; the end-members d — > oo and d = 1 
respectively give BE and FD statistics [HJ [SU] . A dif- 
ferent statistic by Haldane and Wu [311 [32] emerges from 
a generalized Pauli exclusion principle, involving coupled 
interactions between particles in occupancies {rii}. A 
third approach was given as an ansatz by Acharya and 
Swamy [33], and later justified (approximately) by quan- 
tum symmetry [34] , or by exclusion [33J [35] [35] , particle 
kinetics |37j or fractional entanglement |38j arguments. 

The aims of this work are: (i) to report the cor- 
rect, generalized probability distribution for MB, BE and 
FD statistics containing source terms qi, and thereby 
to examine the validity of (0-(|3) and (pSl-fllP); (ii) 
to examine these as special cases of the Polya statis- 
tic, involving "neither independent nor identically dis- 
tributed" ("ninid") sampling [35] gO] ED El M E5; 
and (iii) to report the asymptotic and non-asymptotic 
cross-entropy (divergence) functions and minimum cross- 
entropy (MinXEnt) distributions for Polya systems based 
on the MaxProb principle [35] . It is shown that the gen- 
eralized Polya statistic contains, and therefore justifies, 
the Acharya-Swamy |33j intermediate statistic. 



II. DERIVATIONS 

By enumeration of several simple examples: (a) a non- 
degenerate system; (b) a uniformly degenerate system 



1 In the MaxEnt community, qi is termed the prior probability, 
implying a connection to Bayesian inference; the existence and 
form of such a connection is the subject of much debate [e.g. 

ED [26] ED. 
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FIG. 1: Configurations of (a) Maxwell-Boltzmann, (b) Bose- 
Einstein and (c) Fermi-Dirac ball-in-box models, in which all 
N entities are fully allocated (the gi need not be equal). 



(equal gi); (c) a non-uniformly degenerate system; and 
(d) systems with unequal and qi, it is readily verified 
that the Brillouin distributions ([5|-([7]) - and implicitly, 
the above weig hts - are valid only for the case of 

equal source probabilities for each degenerate sublevel, 
Qi = 9i/G. They do not apply to equal source probabil- 
ities for each level, = s , or to more general forms. 
The Kapur distributions (j8|-(lCM do not satisfy normal- 
ization, X]pe{P} IP = F It is therefore necessary to reex- 
amine the probabilistic basis of each statistic. 

Maxwell-Boltzmann statistics: Firstly, consider a sys- 
tem in which N distinguishable balls (entities) are ar- 
ranged amongst s distinguishable boxes (levels) , in which 
each level contains gi distinguishable sublevels, as shown 
in Figure [T^. Each level has the source probability qi 
of being filled (not necessarily equal), with Yli=i Qi = !• 
We consider each realization of the system, containing m 
balls in each level i, Vi; this in turn is a set of configura- 
tions, or distinguishable arrangements of balls in boxes. 
In counting configurations, we distinguish between rear- 
rangements of balls between boxes, but do not (or need 
not) consider rearrangements of balls within each box 
(see Figure [Tk) . No account is made of the order of alio- 
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eating balls to boxes. Of course, the realizations can be 
chosen in many different ways; the occupancies {rij} are 
selected here for their utility in physical systems. 

Clearly, the probability of an individual arrangement 
of rii balls in the ith level is , and so the probability of 
any specified configuration amongst the levels is Yii=i Qi" 1 ■ 
To obtain PmBi this must be multiplied by the number of 
ways in which the balls can be rearranged amongst the 
levels to give the same realization, ^V!/Ili=i n «" The 
result is the non-degenerate multinomial distribution: 
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Another viewpoint is that if the source probability of the 
jth level is q,, that of the jth sublevel of the ith level is 
lij = Qi/di > V? = 1 , . . . , <7i , assuming equiprobable filling of 
sublevels within a level. The probability of a configura- 
tion amongst the sublevels - i.e. of riy balls in the jth sub- 
level of the ith level, V{i, j} - is then fli=i IlfLi (li/9i) n ' J 

= UU(Qi/9i)^ 9jUnii = UU(Qi/9i) ni - This must be 
multiplied by the number of configurations amongst sub- 
levels, N\ rii=i 9i ni / n i-i again giving (11 1. For equiprob- 
able sublevels qij = G~ 1 , qi = gi/G and (111 reduces to 
the Brillouin form Since the total number of config- 
urations amongst sublevels is G N , this is proportional to 
the degenerate weight ([!]). 

An important conclusion is that when considering a 
generalized source probability qi unrelated to the degen- 
eracy <?,;, the latter becomes irrelevant to the probability 
of each configuration. We can in fact put qi = uii/M, 
where m, € N is a fictitious degeneracy controlled by the 
source probability and M — Yli=i m i> whence: 
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Eq. (12 1 has a strong connection to urn models of prob- 
ability distributions [e.g. [33] HU US], as shown in Figure 
[2] in which balls are drawn from an urn containing M 
balls, made up of m, balls of each level (color) i. In each 
sampling event, a ball is drawn from the urn, its color 
recorded, and is returned to the urn in accordance with 
some rule, until a sample of A~ balls has been recorded. 
In the MB case, the balls are drawn with replacement, 
hence the qi remain constant. As shown in Figure [2] 
the balls are indistinguishable except for their color; the 
weight AM/ Y[t=i n i- now arises from the number of ways 
in which a given set of balls can be sampled, i.e. the num- 
ber of ordered sequences of balls which make up each re- 
alization or type. Urn models thus capture many features 
of sets of alphabetic characters cij drawn from an alpha- 
bet A = {ai}, considered in information theory |4"T1 148j . 
The urn model insight is of benefit to later analyses. 

The above distribution can be generalized to the case 
of unequal sources qij for each sublevel, to give: 
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FIG. 2: An urn model, in which balls are drawn, recorded 
and returned to an urn, until N balls have been sampled. 



with normalization Yli=i Sj=i Qij ~ 1- It is unclear why 
such a complicated arrangement might be needed, since 
each sublevel could be handled as a separate level; but it 
does give a recipe for "coarse graining" if this is desired. 

Bose-Einstein statistics: We now consider the same 
balls-in-boxes problem but with indistinguishable balls, 
as shown in Figure [T]d. The standard argument is that 
there are ( 9i+ ," i _ ) distinguishable ways of arranging m 
indistinguishable balls in gi boxes, Vi, hence the num- 
ber of configurations amongst sublevels in a realization is 
given by the BE weight ^ . For equally probable configu- 
rations, the governing probability is obtained by dividing 
by the total number of configurations amongst sublevels, 
( G+ jv _1 )' g ivin g the Brillouin form 

We now employ the urn model analogy, to consider the 
generalized BE distribution: 
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(12) Using qi — rrii/M and parameter (3 — N/M < 1 gives: 
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with qiN j ' f3~ G N. This is a variant of the multivariate 
negative hypergeometric or Dirichlet compound multino- 
mial distribution [JU |SD], as is the Brillouin BE dis- 
tribution ^ itself. This arises from a complicated sce- 
nario of sampling without replacement, in which two sets 
of balls are drawn from an urn; the conditional distribu- 
tion of the second sample, given the first, is the one of 
interest [42] [50] . For equiprobable sublevels, = gi/G; 
taking = gi and M = G, (14) reduces to the Brillouin 
form and thus is proportional to the BE weig ht (|}. 
For equiprobable levels qi = s^ 1 , thus taking rrii = 1 
and M = s, it yields the non-degenerate BE distribu- 
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tion ( N ~ ) ■ Outside these two cases, (15 1 gives a 
generalized BE distribution with unequal source proba- 
bilities qi , which satisfies normalization and several other 
requirements of a governing probability distribution (see 
later sections). 

Fermi-Dirac statistics: Here we also have indistin- 
guishable balls, but each sublevel contains a maximum 
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of one ball, as shown in Figure ^flp. By standard rea- 
soning, there are distinguishable ways of arranging 
iii £ {0, 1} indistinguishable balls in gi boxes, V«; thus the 
number of configurations amongst sublevels in a realiza- 
tion is given by the FD weight Q . For equiprobable con- 
figurations, the total number of configurations amongst 
sublevels is (^) , giving the Brillouin distribution ^ . 

Applying the urn model analogy gives the generalized 
FD distribution: 
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with the condition m; > n*. In terms of qi and parameter 
P+ = N/M < 1: 



divisor in ( 19 ) arises from the fact that the MB weight is 
normalized differently to the BE and FD weights. Since 
probabilities are always normalized consistently, no such 
discrepancy appears in the Brillouin distributions |[5j|- 
(7L nor indeed in the generalized forms (111, (fl5|) and 



(17 1 derived herein. In consequence, the "corrected Bolt- 
zon" statistic has no independent physical meaning, be- 
ing merely an artefact of the use of weights; indeed, ex- 
cept in the asymptotic limit N — > oo, it gives absurd re- 
sults. This demonstrates an important point, that prob- 
abilistic inference using the MaxProb principle should 
always be based directly on true probabilities P; any con- 
trary method (e.g., one based on weights W) can lead to 
logical inconsistencies [18] , 
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with q l N//3+ £ N and q l N//3 + > m. This and the Bril- 
louin FD distribution Q are variants of the multivariate 
hypergeometric distribution [42 49, 50J, which arises from 
a simple case of sampling without replacement [42 49 , 50J. 
For equiprobable sublevels, this function also reduces 
to the Brillouin form Q, usually expressed as the FD 
weight ([3]), whilst for equiprobable levels, it gives the 
non-degenerate FD distribution (^J . Aside from these 
cases, ( |17[ ) gives a generalized FD distribution with un- 
equal source probabilities qi. 



III. DISCUSSION 

For a, b £ N, the following two inequalities are well 
known: 



«! < b ^ (a + b-l)\ 
(0-6)! ~ " (o-l)! ' 



(18) 



with equality for o> t. Proofs are given in Appendix B. 
These and the weights ([l}-(|3]) yield the inequalities [e.g. 



J FD < < W BE 



(19) 



with equality as gi 3> rii, Vi. If we consider the Brillouin 
distributions (Js])— (JtJ) , one obtains: 



V,D|G < ^MB\G % ^BE\G 



(20) 



The variable direction arises from the terms which con- 
verge to Oi=i 9?'/G N , for which opposite examples can 
be found (see Appendix B); equality in (20) occurs when 
G > N and gi > ra i; Vi. 



Much attention has been paid to (19 1, especially to 



the term Wmb/NI, which has even been labelled a sep- 
arate statistic for "corrected Boltzon" particles |5T] . 
However, comparing ( 19 1-( 20 ) , it is evident that the Nl 



IV. THE POLYA DISTRIBUTION AND 
MAXPROB 

We now consider the Polya urn model (Figure [2]) , in 
which a ball is sampled from an urn, recorded, and re- 
turned to the urn; in addition, c £ Z balls of the same 
color (level) are also added to the urn. The probability of 
drawing the realization {n^} from the urn after N draws 
is [391 Sol MM H3H3H: 
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(21) 

The Polya urn thus describes a simple, analytic case 
of "neither independent nor identically distributed" 
( "ninid" ) sampling, of which the bivariate form is widely 
used in studies of biological contagion [e.g. |4"3"]. 

The Polya model (21 1 can be written in the following 
forms, valid respectively for c = 0, c > and c < 0: 
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These distributions encompass the MB, BE and FD gov- 
erning distributio ns (fl2j ), (14 1 and ( [16]). Simplified forms 
of @ and fl23])-(|24fl were given bvp6"ll37]. In light of 
the Polya model, we can set aside the complicated urn 
model for the multivariate negative hypergeometric dis- 
tribution (the generalized BE statistic (14l), to offer the 
following explanations for each statistic: 

• MB statistic: Polya urn with c=0 (sampling with re- 
placement); 

• BE statistic: Polya urn with c=l (sampling with dou- 
ble replacement); 
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• FD statistic: Polya urn with c=-l (sampling without 
replacement) ; 

It is curious that the last two representations are equiv- 
alent to the allocation of indistinguishable balls to dis- 
tinguishable boxes (with a further restriction in the FD 
case). Whilst the full implications of this insight are not 
clear, it does provide an alternative philosophical foun- 
dation for the construction of quantum statistics. 



The Polya distributions ( 23 )-( 24 1 can be further ana- 
lyzed using an extension of Boltzmann's principle (HI: 
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where D is the generalized cross-entropy (directed diver- 
gence or negative relative entropy) function defined by 
the MaxProb principle. Applying ( p5| in the asymptotic 
limit N — ► oo to ([22t-(p4|, using pi = m/N, q t = mjM 
and (3 = N/M, gives: 
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(28) 

with restriction |c| < qi//3pi,\/i (whence |c| < mj/ni,Vi) 



for c < 0. Note (26 1 is the Kullback-Lieblcr cross-entropy 



[52] . Proof of these results based on Sanov's theorem is 
given elsewhere |45j : a more simplistic calculation may 
be conducted by applying Stirling's approximation to 
all factorials and gamma functions. Both forms of the 
Polya cross-entropy (27 1- (28) have been shown to satisfy 



many important information-theoretic properties, includ- 
ing non-negativity, lower semicontinuity, convexity, par- 
tition inequality and the conditional limit theorem [45 . 
Applying the MaxProb principle, the relevant cross- 



entropy (26 1, (27 1 or (28 1 can be minimized (MinXEnt) 



subject to the following constraints: 



j> = l. ( 29 ) 
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J2Pif" = (fr), r=l,...,R, (30) 



where f r i is the value of the function f r in the ith level 
and (f r ) is its mathematical expectation. Expremiza- 
tion by the Lagrangian method gives the inferred "most 
probable" distribution for each statistic: 
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where A r ,r = 0,...,R are Lagrangian multipliers, with 
constants absorbed into the Ao terms. For equiprob- 
able sublcvcls qi = gi/G, sample size (3 — N/G and 
the above values of c, the cross-entropy functions p6| - 
(28 1 and distributions (31 1- ( 33 > reduce to those given 
by Brillouin [51 [5]; these in turn simplify - up to con- 
stant factors - to the commonly used entropy functions 
and distributions obtained from the weig hts 2 0-j3) 

[DISIEISIEIISIISIIIDHIIIISI. 

For N <C cxd, the non-asymptotic cross-entropy func- 
tions and distributions for the above systems can also be 
derived using Boltzmann's principle (251, as reported in 



Appendix A. These encompass the non-asymptotic MB, 
BE and FD entropy functions and distributions reported 
previously |T6j [17] . 

To summarize, generalized MB, BE and FD statistics, 
containing unequal source probabilities for each level, can 
be obtained using the combinatorial definition of cross- 
entropy ( 25 1 in conjunction with probability distributions 
generated using urn models. All such cases emerge as 
special instances of the Polya urn model, involving a sim- 
ple scenario of urn modification during sampling. Of the 
systems examined, only the MB system in the asymp- 
totic limit yields the Kullback-Leibler cross-entropy ( 26 1 



|52j : in consequence, the latter is not universal in appli- 
cation, and cannot be used to infer the most probable 
realization of BE, FD or Polya systems except in special 
limiting cases (nor even of the MB system away from 
the asymptotic limit N — > oo). This conclusion con- 
tradicts the dominant viewpoint in information theory, 
in which the Shannon entropy |47j and Kullback-Liebler 
cross-entropy [55] are considered to be the paramount 
functions for inductive reasoning ( "statistical inference" ) , 
based on their axiomatic and information-theoretic defi- 
nitions [c.f. 1531 154"] . This conclusion is, however, a neces- 
sary consequence of Boltzmann's principle. 



Eqs. (31 )-( 33 ) can be amalgamated as: 
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(34) 



2 Physicists typically consider constraints on normalization | |29| l 
and energy fu = Ei, (/i) = (E) in ( |30| , giving multipliers 
Ai = 1/kT and Ao = —fi/kT, where k is the Boltzmann constant, 
T is absolute temperature and fi can be viewed as a chemical 
potential. 
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This can be compared to the ansatz presented by Acharya 
and Swamy [33] . in the present notation: 



Pasa — 



9i/N 



(35) 



in which a = 0, 1 and —1 respectively give MB, BE and 
FD statistics, and — 1 < a < 1 their intermediate interpo- 
lation. The Polya representation (34 1 thus contains the 



Acharya-Swamy intermediate statistic ( 35 ) as a special 
case. It therefore justifies this statistic using a "ninid" 
urn sampling scheme, an explanation rather different to 
that of previous studies [3 ^ 13 ^ 1551, |3"6l |3"?1 155] . 



V. CONCLUSIONS 

Generalized probability distributions for MB, BE and 
FD statistics, with unequal source probabilities qi for 
each level i, are derived by extension of the distribu- 
tions given by Brillouin [SJ [5] , using an urn model anal- 
ogy. These are shown to be special instances of the Polya 
urn model, involving a simple scenario of urn modifica- 
tion during sampling ("ninid" sampling). The resulting 
cross-entropy functions and MinXEnt distributions are 
derived by Boltzmann's principle (the MaxProb method). 
The general form of the Polya distribution contains the 
Acharya-Swamy |33j intermediate statistic, as an exact 
result without approximation. 

We also show that the "corrected Boltzon" statistic 
based on Wms/JV! has no independent meaning, but is 
an artefact of the use of weights instead of probabili- 
ties in Boltzmann's principle. Furthermore, since the 
Polya cross-entropy function in general differs from the 
Kullback-Liebler cross-entropy [52], the latter does not 
give the most probable realization of BE, FD or Polya 
systems except in special limiting cases. 

As a final comment, the statistic of Haldane and 
Wu [31] [32] considers intermediacy as a consequence of 
changes in occupancies {n^} due to interactions. The 
Acharya-Swamy-Polya model considers intermediacy due 
to changes in the source probabilities {qi} - hence in the 
degeneracies {gi} - during sampling. Either approach is 
justifiable on physical grounds; their common features 
and possible synthesis deserve further investigation. 
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VI. APPENDIX A: NON-ASYMPTOTIC 
DISTRIBUTIONS 



The non-asymptotic forms of the Polya cross-entropy 
functions can be obtained directly from Boltzmann's 



principle ( 25 1 , without taking the asymptotic limit N 



oo. Such statistics are required for systems of small num- 
bers of non-interacting entities 16, 17, 18, 19 . From ( 25 1 
and the governing distributions (11), (15) and (17 1, the 
functions are: 



c=0 



N 



i=i 



-D 



(N) = 
Polya iu 
c>0 JV 



l^j M T(N/0c) 
r(N/pc + N) 

T(q t N/f3c + Pl N) 



In 



( Pl N)\ T( qi N/pc) 



(36) 
(37) 



-D 



(N) _ 
Polya 



c<0 



i=l 
hi 



Pi In 



N\ T(l - N/Pc - N) 



r(l - N/(3c) 
r(l - qi Nl(3c) 



(p x Ny.T(l- qi N/Pc- Pl N) 



(38) 



in which constant terms are brought inside the sum using 
(JggT). Applying MinXEnt subject to p9|-p0| gives the 



inferred "most probable" distribution for each statistic: 



PPolya,', 
c=0 



Ia«i-Ar°-X)A r / ri } 
r=l J 



(39) 



P Polya,, = T7 A ^ 
c>0 JV 



N\ T(N/f3c) 



N T(N/f3c + N) 

R 

A( qi N/[3c + p* olyai N~ 1)-A§> -J2*rfr 



Ppoly 
c<0 



c>0 r —i 



1 ! f 1 N\T(l-N/(3c-N) 
N \N r(l - N/0c) 



(40) 



A(- qi N/Pc- P * olyai N)-\ c <° 



c<0 



R . 

^ ^rfri r 
r=l ' 



(41) 



where X r ,r = 0, ...,R again are Lagrangian multipliers, 
A(x) = ip(x + 1) = y is a shifted digamma function, de- 
fined for convenience, and A _1 (y) = -0 _1 (y) — 1 is its 
inverse. None of these functions allow explicit factoriza- 
tion of a partition function Z = e A ° ; the last two forms 
are also implicit in pf . 
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VII. APPENDIX B: INEQUALITIES 



For b > the terms in ( 18 1 can be written: 

a(a-l)...(a-b+l) (42) 
(a + b-l)(a + b-2)...a (43) 



(a -b)\ 
(a + b-l)\ 



(a-l)\ 



As a/b — > oo, each factor in (42l-(43l will approach 
respectively from below or above 
verge to a b in accordance with (18 1; (19 1 follows 
From (El-pl, both 



Both forms thus con- 
□ 



be\g and ¥ FD \ G converge 
Vmb\g- However, the sign of the inequalities can vary. 
This is illustrated by the set of realizations for the Bril- 
louin case ./V = 4, s = 3 and gi = 3,Vi, for which the 
probabilities ^ are listed in Table [I] (some realiza- 
tions are excluded from the FD set). Depending on the 
realizations, the order can be FprjiG < ^mb\g < ~ 



to 



FD\G 



BE\G, 



™FD\G > ^MB\G > ^BE\G or ™FD\G < ^MB\G > ^BE\G- 

In this and many other examples calculated, we did not 
observe Pfd|g > ^mb\g < ^be\g, although this result 
is not proven. 



Realization 
[ni,n 2 ,n 3 ] 




AIIV 




AIIV 




[1, 


1, 


2] 


0.21429 


> 


0.14815 


> 


0.10909 


[1- 


2, 


1] 


0.21429 


> 


0.14815 


> 


0.10909 


[2, 


1, 


1] 


0.21429 


> 


0.14815 


> 


0.10909 


[2- 


2, 


0] 


0.07143 


< 


0.07407 


> 


0.07273 


[2- 


0, 


2] 


0.07143 


< 


0.07407 


> 


0.07273 


[o. 


2, 


2] 


0.07143 


< 


0.07407 


> 


0.07273 


[1, 


3, 


0] 


0.02381 


< 


0.04938 


< 


0.06061 


[b 


0, 


3] 


0.02381 


< 


0.04938 


< 


0.06061 


[3- 


1, 


0] 


0.02381 


< 


0.04938 


< 


0.06061 


[3- 


0, 


1] 


0.02381 


< 


0.04938 


< 


0.06061 


[0. 


1, 


3] 


0.02381 


< 


0.04938 


< 


0.06061 


[o. 


3, 


1] 


0.02381 


< 


0.04938 


< 


0.06061 


[4- 


0, 


0] 


NA 




0.01235 


< 


0.03030 


[0, 


4, 


0] 


NA 




0.01235 


< 


0.03030 


[o. 


0, 


4] 


NA 




0.01235 


< 


0.03030 



TABLE I: Probabilities of each realizations for N 
and g = [3,3,3] (NA=not applicable). 
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